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The electronic dispersion for Bi2Sr2CaCu2 08+i has been determined from angle-resolved pho- 
toelectron spectroscopy (ARPES). From this dispersion we calculate the entropy and superfluid 
density. Even with no adjustable parameters we obtain an exceptional match with experimental 
data across the entire phase diagram, thus indirectly confirming both the ARPES and thermo- 
dynamic data. The van Hove singularity is crossed in the overdoped region giving a distinctive 
linear-in-T temperature dependence in the superfluid density there. 
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The generic doping dependence of the thermodynamic, 
electrodynamic and transport properties of high- Tc su- 
perconductors remains a puzzle despite many years of 
study. Their unusual behaviour is often taken to be a 
signature of exotic physics yet it should be related to 
the electronic energy-momentum dispersion, obtained for 
example from angle-resolved photoemission spectroscopy 
(ARPES). These studies indicate the presence of an ex- 
tended van Hove saddle-point singularity[l] situated at 
the (0,7r) point together with a normal-state pseudogap[2] 
as common features in the electronic band structure of 
the cuprate superconductors. The pseudogap exhibits a 
reduction in the density of states (DOS) at the Fermi 
level which is believed to develop into a fully nodal gap 
at low temperature [3i|. 

In theories based on the so-called van Hove scenario Q 
the superconducting (SC) transition temperature, Tc, 
is enhanced by the close proximity of a van Hove sin- 
gularity (vHs). These theories assume that the vHs 
sweeps through the Fermi level at around optimal dop- 
ing {p = 0.16), indeed causing the peak in Tc{p). How- 
ever, ARPES measurements of the band structure of Bi- 
22010 show that the vHs crosses the Fermi level in the 
deeply overdoped side of the phase diagram. For a bilayer 
cuprate like Bi-2212 the weak coupling between the lay- 
ers splits the bands near the (tt, 0) points into an upper 
antibonding band and a lower bonding band. ARPES 
measurements performed on Bi-22120 suggest that the 
antibonding band vHs crosses at around p = 0.225 where 
Tc ~ 60K i.e. near the limit for overdoping in this mate- 
rial. The vHs crossing should profoundly affect all phys- 
ical properties. 

In this work we have used the thermodynamic proper- 
ties as a window on the electronic structure to indepen- 
dently check the main ARPES results. Using an ARPES- 
derived energy dispersion we have calculated the doping 
and temperature dependence of the entropy and super- 
fluid density of Bi-2212. All details for our calculations 
are taken directly, and only, from ARPES in order to de- 
termine the implications of this data. Our calculations 



confirm the ARPES results, giving a consistent picture 
of the thermodynamic and electrodynamic properties in 
terms of a proximate vHs. 

We assume Fermi-liquid-like, mean-field, weak- 
coupling physics in spite of indications, or expectations, 
to the contrary. In defense of our approach (i) the ther- 
modynamics at low T is dominated by the nodal regions 
of the Fermi surface where quasiparticles are well defined 
and long lived; and (ii) the Wilson ratio relating spin 
susceptibility to S/T (where S is the electronic entropy) 
is almost exactly that for nearly free electrons across a 
wide range of doping and temperature?]. Moreover, the 
non-mean-field BCS-like behaviour is generally inferred 
from the unusual 2A/kBTc ratio which grows with un- 
derdoping. But it has always been our view[8], and is 
now confirmed!^, [loi] that the large and growing energy 
gap used here is the (7r,0) pseudogap, not the SC gap. 
Once the pseudogap is properly included in the problem 
then 2A/kBTc is well behaved. We do not treat fluctua- 
tions which are confined to Tc ± 15K[llj and are a minor 
embellishment. 

For Bi-2212 we employ a 2D bilayer dispersion ek pro- 
vided by the authors of ref . @ , which was obtained from 
tight binding fits to high-resolution ARPES data. The 
DOS per spin at energy E is given by 

k 

The entro py p er mole S for weakly interacting fermions 
is given bv[12l| 

S=-2Rj [/ In / + (1 - /) In (1 - f)]N{E)dE (2) 

where / is the Fermi-Dirac distribution function and R is 
the gas constant. The chemical potential /^(T) is calcu- 
lated self-consistently such that the carrier concentration 
n is T-independent. n is given by: 

n = {2/Va) I f{E)N{E)dE (3) 
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FIG. 1: (Color) The DOS calculated from the Bi-2212 bilayer 
dispersion determined by ARPES measurements (black). 
Also shown is a 20meV pseudogap at OK (green) and lOOK 
(blue), and a 20meV SC gap at OK (red). Inset: the Fermi 
surface in the {kx,ky) plane showing the angle 9. 



Where Va is the atomic volume per formula unit. 

The Fermi surface (FS) in the 1st Brillouin Zone is 
shown in the inset to Fig. [TJ The pseudogap first forms 
on the FS near (tt, 0) leaving ungapped Fermi arcsfisj 
between. With decreasing temperature the Fermi arcs 
narrow such that the gap seems to become nodal at T = 
0. We therefore adopt a pseudogap of the form 



where 



Eg^fnax COS 



{0 < Oo) 



> 



%) (4) 



otherwise 



TT ( , I T 

70 = — I 1 — tanh 



(5) 



and T* = Eg^max/ks- is the angle shown in Fig.[TJ 

Eqn. [5] models the observed temperature dependence 
of the Fermi arc lengthQ. At T = 0, 6^0 = 7r/4 and the 
pseudogap is fully nodal. As T rises, decreases result- 
ing in a 'filling-in' of the pseudogap and the growth of the 
Fermi-arcs. This model is based on results by Kanigel et 
aZ. Q that show the Fermi-arcs collapsing linearly as a 
function of T/T* , extrapolating to zero as T ^ 0. How- 
ever we note an important feature of our model. The 
Kanigel data shows the pseudogap opening abruptly at 
T = T* . A pseudogap which fills completely at T* would 
result in a jump in the specific heat coefficient 7 at T*, 
which is not observed. The smooth evolution of the tanh 
function in Eqn. [5] overcomes this problem. The pseudo- 
gap is states-non-conserving i.e. unlike the SC gap there 



is no pile up of states outside the gap (see Fig. [T]). This is 
implemented by eliminating states with energies E < Eg 
from the summations. 

Fig. [T] shows the DOS calculated from the bilayer dis- 
persion. The bonding and antibonding band vHs's are 
clearly visible with the former 105meV below the latter. 
Also shown is a 20meV pseudogap at T = OK and lOOK 
illustrating the gap filling with temperature. The gap 
node is pinned to the chemical potential at all T. 

The entropy in the SC state has been modelled using a 
c?-wave gap of the form Ak = ^AqQw where = cos — 
cos ky. The dispersion in the presence of the SC gap is 
given by i?k = v^k+^k ^^'^ ^0 (T) is determined from 
the self-consistent weak-coupling BCS gap equation [lij 
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We adopt a pairing potential of the form Vkk' — ^5k5k' • 
The amplitude, V , is assumed to be constant (=125meV) 
up to an energy cut-off, Wc, chosen such that matches 
the experimentally observed value. The pseudogap is not 
included in the process of calculating Ao(T). 

The superfluid density, ps, is proportional to the in- 
verse square of the penetration depth given bv[l^ 
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The summation in Eqs. [5] and [7] is performed over both 
the bonding and anti-bonding bands and A is assumed 
to be the same for both bands [isj. 

The data points in Fig.[2Ua) show the normal- and SC- 
state entropy data of Loram et aZ0. Fig.[2tb) shows the 
absolute entropy calculated, as described, from the dis- 
persion with no fitting parameters. We have merely spec- 
ified the location of Ep relative to E^hs at two points 
only: in the overdoped region from Kaminski et aj.[6| 
and in the underdoped region from Kordyuk et aL[l6l| 
and interpolated between. The overall T- and doping- 
dependence of the experimental data is reproduced su- 
perbly, with absolute values just a factor of 3/4 lower. 
This is our first main result. S/T rises with doping and 
reaches its maximum at the vHs (as observed also in 
La2-a;Sra;Cu04[3|)- A similar effect is seen in the spin 
susceptibility and Knight shift Tl*!. In Fig. [2{a) we 
have rescaled the computed entropy by the constant fac- 
tor of 4/3 and refined the fit by using Ep and Eg^max 
as fitting parameters. These refinements do not alter the 
overall behaviour and are tightly constrained. For exam- 
ple, the normal-state fits to the four most overdoped data 
sets have been obtained by adjusting a single parameter, 
namely Ep, as are the high-T asymptotes for all data 
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FIG. 2: (Color) (a) Refined normal-state (black) and SC-state 
(red) fits to the electronic entropy data of Loram et al^ for 
Bi-2212. For clarity every 20th data point only is shown. 
Each curve represents a different doping level from p = 0.129 
to 0.209. (b) Unrefined absolute entropy curves with no fitting 
parameter. 
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FIG. 3; (Color) (a) Ep - E^hs (up triangles). The vHs 
will cross Ef at p ~ 0.22. Pairing potential energy cut- 
off ujc (down triangles). Pseudogap magnitude Eg (squares) 
and a fit given by Eqn. [8] The measured Tc multiplied by 
2.14fcs(circles). (b) Maximum gap as measured from the DOS 
at lOK (circles). SC gap measured from the DOS at lOK in 
the absence of a pseudogap (diamonds). Calculated super- 
fiuid density at lOK (up triangles). Measured low-T Bi-2212 
superfiuid density [2l| (down triangles). 



sets. The Fermi level in the most overdoped fit is only 
8meV above the antibonding band vHs. 

As the doping decreases the vHs recedes from Ep re- 
sulting in a decrease in the number of states within fc^T 
of /i and a corresponding reduction in entropy. How- 
ever as the doping is further reduced the recession of 
the vHs from Ep is no longer able to account for the 
observed decrease in entropy alone and the second ad- 
justable parameter, the pseudogap magnitude Eg, is in- 
troduced. This results in the progressive downturn in the 
normal state S/T as temperature decreases. The deduced 
values of Ep ~ E^hs and Eg are plotted versus doping 
in Fig. Efa) along with the measured Tc- The doping 
level has been determined from the empirical relation 17] 
p = 0.16 ±0.11^/1 - Tc/Tc^max- The fits suggest that the 
antibonding vHs will cross Ep near p = 0.22 in full agree- 
ment with recent ARPES studies 6] on Bi-2212 where the 
crossing occurs at p = 0.225. This is our next key result. 

The pseudogap is observed to open at critical do pin g 
Pcrit = 0.188 in agreement with previous analvses[^ Il8|. 
Eg has been fitted with the following equation 



T*{p) = Eg/ks = Tq (1 - p/pcrit)^~ 



(8) 



with Tq = 443. 7K and a — 0.317. These values agree 
with the results of Naqib et al. [l^ who determined T* (p) 
of YBCO from transport studies. A fit to their data gives 
T* = 510K and a = 0.2. The sublinear behavior of T*(p) 
is expected if Pcrit is a quantum critical point [20| . 

From the SC state fits the energy cutoff, ojc, is, sur- 
prisingly, found to be linearly related to Ep — E^hs- In 
particular Wc(?Tiey) = 10.96+0.637{Ep — EyHs) with cor- 
relation coefficient R = 0.99945. The rapid fall in ujc is 
suggestive of magnetic or magnetically enhanced pairing, 
rather than phononic. 

Fig. [S^b) shows the magnitude of the combined SC 
gap and pseudogap, Amax, measured from the calculated 
DOS at lOK. Afnax increases with decreasing doping just 
as observed from ARPES [l^l, tunnclling[23] and Raman 
scattering [2J]. Also plotted is the SC gap magnitude A, 
determined by setting Eg = and measuring the gap in 
the calculated DOS at lOK. The magnitude is smaller 
than typically observed because of the weak coupling as- 
sumption for which 2A/kBTc = 4.28. The gap. A, rises 
and falls in conjunction with the observed Tc- Note that 
the experimentally observed monotonic increase in the 



gap magnitude with decreasing doping is here seen to 
be associated with the pseudogap, and not the SC gap as 
generally believed. The behaviour here is consistent with 
the two-gap picture presented by Deutscher[25'| and more 
recently by Le Tacon et al. llCl| but has been a feature of 
our work for a long timefsl-jisl. [26j. 

Using the parameters obtained from the entropy fits 
the superfluid density has been calculated using Eqn. 
[7] with no further adjustable parameters, and is shown 
in Fig. m^a). For comparison (and in the absence of 
data for Bi-2212) Fig. [l^b) shows the a6-plane super- 
fluid density of La2-xSr.T;Cu04 (La-214) determined by 
ac-susceptibility measurements [271 grain-aligned sam- 
ples. There is excellent agreement. The increasing lin- 
earity of Ps{T) with overdoping can now be understood 
in terms of the approach to the vHs where full linearity 
occurs. (The crossing of the vHs in La-214 can also be in- 
ferred from the maximum in the entropy at p = 0.24(3].) 
The opening of the pseudogap leads to the strong re- 
duction in ps observed below p — 0.19. This is clearly 
illustrated by the plot of Ps(lOK) vs p in Fig.Efb). The 
overall doping dependence and absolute magnitude of 
Ps(lOK) concurs almost exactly with experimental data 
for Bi-2212[2lj also shown in Fig.[3tb). This is our third 
key result. We recall that no fitting parameters are used 
in Eqn. [71 It is remarkable that S/T and Ps{T) are so sim- 
ilar in La-214, Bi-2212 and indeed Yi_2;Caa;Ba2Cu307_5, 
despite the significant differences in bare band structure. 
The renormalised dispersion near Ep seems to lead to a 
universal phenomenology which calls for theoretical ex- 
planation. 

In the underdoped data the downturn seen at low T 
and p in the calculated ps (T) curves arises from the clos- 
ing of the Fermi arcs and is not observed in the exper- 
imental data which show an upturn at low T and p. A 
similar downturn would occur in the condensation en- 
ergy, which again is not observed. To us this indicates 
that the Fermi arc picture is, at best, incomplete. We 
will discuss this elsewhere. 

In summary, we have calculated the entropy and su- 
perfluid density of Bi-2212 directly from an ARPES- 
derived energy-momentum dispersion. The temperature 
and doping dependence of both the entropy and super- 
fluid density can be fully explained by the combined ef- 
fects of a proximate vHs and the opening of a normal- 
state pseudogap. These results provide indirect confirma- 
tion of both the thermodynamic and low-energy ARPES 
data. Fits to Bi-2212 entropy data indicate that the 
antibonding vHs crosses the Fermi level near p = 0.22 
in agreement with recent ARPES results. The super- 
fluid density calculated using no adjustable parameters 
shows excellent agreement with experimental data and 
exhibits a distinctive overall linear-in-T behaviour at the 
vHs. The universal renormalised phenomenology in the 
various cuprates, despite their differences in bare band 
structure, is a key conclusion that demands theoretical 
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FIG. 4: (Color) (a) Superfluid density computed using para- 
maters from the entropy fits in Fig. [2] The arrow points in 
the direction of increasing doping. Curves/data in which the 
pseudogap is present are shown in blue, (b) Superfluid den- 
sity of La-214 obtained from ac-susceptibility measurements 



explanation. It also remains a theoretical challenge to 
understand why the Fermi-liquid approach is so success- 
ful in a strongly-correlated system. 

Thanks are due to A. Kaminski for the dispersion pa- 
rameters and also Prof N.W. Ashcroft for much discus- 
sion and many helpful comments on this work. 
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